Response of the Strongly-Driven Jaynes— Cummings Oscillator 
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We analyze the Jaynes-Cummings model of quantum optics, in the strong-dispersive regime. In 
the bad cavity limit and on timescales short compared to the atomic coherence time, the dynamics 
are those of a nonlinear oscillator. A steady-state non-perturbative semiclassical analysis exhibits a 
finite region of bistability delimited by a pair of critical points, unlike the usual dispersive bistability 
from a Kerr nonlinearity. This analysis explains our quantum trajectory simulations that show 
qualitative agreement with recent experiments from the field of circuit quantum electrodynamics. 

PACS numbers: 42.65.Pc, 42.50.Dv, 42.50.Pq, 03.65.Sq 
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The Jaynes-Cummings (JC) Hamiltonian provides a 
quantum model for a two-level system (qubit) interact- 
ing with a quantized electromagnetic mode. It is widely 
applicable to experiments with natural atoms ^|H] as 
well as for solid-state 'artificial atoms' [9l411j. The JC 
Hamiltonian can be diagonalized analytically, but in the 
presence of a drive ^(i) and dissipation the open-system 
model becomes non-trivial, with the effective behavior 
depending strongly on the specific parameter regime. 
The case where the cavity relaxation rate n greatly ex- 
ceeds the two-level dissipation and dephasing rates 7, 70 
is known as the bad cavity limit. The strong dispersive 
regime [T^ [13] of the JC model describes the situation 
that the presence of the qubit causes the cavity frequency 
to be shifted by an amount x much greater than the cav- 
ity lincwidth. Recent experiments [14 , which operate 
in both the strong dispersive regime and the bad cav- 
ity limit, show a nontrivial response under conditions 
of strong drive, arousing interest due to its usefulness 
for high-fidelity qubit readout. A characteristic feature 
of the JC model is that for very high excitation num- 
ber A^ ^ 1, the excitation number-dependent shift obeys 
xi^) ~^ 0: the transition frequency returns to the bare 
cavity frequency. In the presence of dissipation this hap- 
pens effectively when x(^) ^ '^i a-nd for all larger N the 
response of the system is linear with respect to the drive. 
We describe this behavior as setting in at an excitation 
number A'barc; with the definition x(-^barc) ^ n. In the 
strong dispersive regime we have A'barc S> A^ait, where 
A^crit as usual denotes the excitation level where the dis- 
persive approximation breaks down (defined below). The 
latter inequality has an important consequence for the 
theory: a perturbative expansion in the small parame- 
ter N/Ncrit, typically useful [15] in the dispersive regime, 
is not applicable for the interesting regime A^ > A^baro 
where the system regains the linear response. 

In this paper we consider the JC model under very 
strong driving, such that A^ :s> NcHt- Our main result 
is that there exists a threshold drive fc2 at which the 
photon occupation increases by several orders of mag- 
nitude over a small range of the drive amplitude. We 
perform both quantum trajectory simulations and a non- 



perturbative semiclassical analysis, including the drive 
and the cavity damping. Our results are in qualita- 
tive agreement with recent experiments [14] for a circuit 
quantum electrodynamics (QED) device |16j containing 
4 transmon ^7\ [18] qubits, demonstrating that the JC 
model captures the essential physics despite making an 
enormous simplification of the full system Hamiltonian. 

The behavior of the JC nonlinearity goes beyond the 
Kerr nonlinearity that is often considered. Dispersive 
bistability [19] from a Kerr nonlinearity is well-known 
in atomic cavity QED |20j. It has been implemented 
in the solid state via the nonlinearity of a Josephson 
junction |21j . and in the circuit QED architecture has 
produced high-fidelity readout of qubits [HJ [53] . Simi- 
lar schemes have been implemented with nonlinear mi- 
cromechanical resonators |;24]. However, unlike the Kerr 
anharmonicity, the JC anharmonicity does not remain 
constant but rather diminishes toward zero as the cav- 
ity occupation is increased. As a result, for sufficiently 
strong drive the response of the JC model must return to 
the linear response of the bare cavity. Instead of coher- 
ent driving, an alternative way to saturate the qubit and 
cause the JC system response to return to the bare cav- 
ity response is to couple the system to a bath at elevated 
temperature, as has been investigated theoretically [5S] 
and experimentally |26l . 

Our analysis is applicable to any experiment that can 
reach the strong dispersive limit and drive sufficiently 
strongly, but for concreteness we adopt the language and 
typical parameters of the field of circuit QED. We write 
the driven JC Hamiltonian 



H = Wcfl'a + -^<yz + 5(a<7+ + a cr_) 
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(a -I- a)) cos{u}dt), 
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with cavity frequency cJc/Stt, qubit frequency a;q/27r, cou- 
pling strength g, drive amplitude ^ and drive frequency 
Wd • Operating in the strong-dispersive bad-cavity regime 
defines a hierarchy of scales 



7, 70 < '^ < 9^/S < (/ < (5 < Wc, 
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Figure 1. Transmitted heterodyne amplitude \{a)\ as a func- 
tion of drive detuning (normalized by the dispersive shift 
X ~ g'^ /5) and drive amplitude (normalized by the ampli- 
tude to put n = 1 photon in the cavity in linear response, 
^1 = fi:/\/2). Dark colors indicate larger amplitudes, (a) Ex- 
perimental data ^22;, for a sample with a cavity at 9.07 GHz 
and 4 transmon qubits at 7.0, 7.5, 8.0, 12.3 GHz. All qubits are 
initialized in their ground state, and the signal is integrated 
for the first 400 ns ~ 4/k after switching on the drive, (b) Nu- 
merical results for the JC model of Eq. (|8|, with qubit fixed to 
the ground state and effective parameters 5/2-k — —1.0 GHz, 
g/2n = 0.2 GHz, K/27r = 0.001 GHz. These parameters are 
only intended as representative numbers for circuit QED and 
were not optimized against the data of panel (a). Hilbert 
space is truncated at 10,000 excitations (some truncation arti- 
facts are visible for the strongest drive) , and results are shown 
for time t — 2.5/k. 



where S = uiq — uJc is the qubit-cavity detuning. We can 
make the standard transformation (28) to decouple the 
qubit and cavity 
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dropping terms from the transformed drive that are sup- 
pressed as ©(n"^/^) and 0{g/5). The resulting Hamil- 
tonian would be trivial were it not for the fact that the 
transformation T, defined by 



T = exp[-6i(4iV)"i/2(acr+ + aV_)], 
sin e = -257V1/VA, cos 9 = (5/A, 

depends on N , the total number of excitations: 



N = a^a+-^ + 



(4) 
(5) 
(6) 



(7) 



For photon decay at rate k we can write the decoupled 
quantum master equation after dropping small terms. 



p = -i[H, p] + K{[ap, flT] + [a, paT])/2, 



(8) 



which we integrate numerically in a truncated Hilbert 
space using the method of quantum trajectories, after 
making the rotating wave approximation (RWA) with 



respect to the drive. The experiments we wish to de- 
scribe are performed on a timescale short compared to 
the qubit decoherence times 7~^, 7!^ and we therefore 
treat cr^ as a constant of motion. The remaining de- 
gree of freedom constitutes a Jaynes-Cummings oscilla- 
tor. Note that the qubit relaxation and dephasing terms 
that we have dropped involve the a± and CTz operators 
and would transform in a nontrivial way under the de- 
coupling transformation T [TS]. The results of the nu- 
merical integration for cr^ = — 1 are compared with re- 
cent experimental data [57| in Fig. Ill where we show 
the average heterodyne amplitude \{a)\ as a function of 
drive frequency and amplitude. Despite the presence of 
4 qubits in the device, the fact that extensions beyond a 
two-level model would seem necessary since higher levels 
of the transmons are certainly occupied for such strong 
driving |29j, and despite the fact that the Rabi Hamilto- 
nian might seem more appropriate for such large photon 
occupation, vN g ~ Wc, nevertheless the JC model quali- 
tatively reproduces the features of the experiment ,30j . In 
particular, for weak driving we see a response as expected 
at the dispersively shifted cavity frequency Wc — Xj with 
X = g^ 1^-, which shifts towards lower frequencies as the 
drive increases. For stronger driving a dip appears in the 
response, which we interpret as a consequence of plotting 
the ensemble-averaged | (a) | in the classically bistable re- 
gion, as we discuss below. For increasing drive the dip 
shifts to lower frequencies; finally for the strongest driv- 
ing, the response becomes centered at the bare cavity fre- 
quency a;c/27r and is single-peaked and extremely strong. 
We note that both the experiment and numerical inte- 
gration are terminated at a transient time of only a few 
cavity lifetimes, and we have checked that the numerical 
response is substantially different for the steady state. 

Due to the large number of photons in the system, 
it is possible to form a semiclassical model, similar to 
Refs. [3TJ [35]. This will be a good approximation in 
the case that ratio of the anharmonicity of the dispersive 
Hamiltonian to the decay rate is such that the N ~\ -^ N 
photon peak overlaps well with the iV o A^ -I- 1 photon 
peak, N ^ iVgc, where TVgc = g'^ 1 1^^^ (for the parame- 
ters of Fig. Ill A'sc = 1.6). In the opposite limit we will 
see photon blockade and associated effects, as in Ref. [33] . 
Recently it was shown that it is possible to have a coexis- 
tence of both the semiclassical and quantum solutions for 
a certain range of parameters of the system and drive [34j . 
The semiclassical model will remain valid for N > A^critj 
where a perturbative expansion of the Hamiltonian (|3| in 
terms of N/Nc^t fails to converge, where A'crit — 5 /Ag. 
We rewrite the Hamiltonian Eq. pi using canonical vari- 



ables X = 




and P = i-v 
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(Tz) + ^X COii{uj^t) 



^^2g^X^ + P^+az) + S^. 
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The semiclassical approximation consists of treating X 



and P as numbers, and the effect of cavity relaxation 
is incorporated through a damping term proportional to 
K. We solve for the steady state, treating X'^ + P^ as a 
constant (thus we ignore harmonic generation), giving a 
nonlinear equation for the amplitude A = y/X^ + P^ 



A' 
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(10) 



with amplitude-dependent frequency shift x(^)i given by 



X{A) = <y. 



v/2g2(A2+a,) + <52- 



(11) 



This reproduces for small driving the usual dispersive 
shift x(0) — ±5^/(5 and for strong driving shows the sat- 
uration effect lim^_>.oo x(^) = 0- 

The solution of Eq. ([lO| is plotted in Fig.[2]for the same 
parameters as in Fig. [Ij For weak driving the system re- 
sponse approaches the linear response of the dispersively 
shifted cavity. Above the lower critical amplitude ^ci 
the frequency response bifurcates, and the JC oscillator 
enters a region of bistability. We denote by Ci the point 
at which the bifurcation first appears. Dropping terms 
which are small according to the hierarchy of Eq. (l2]) , this 
point occurs at 



e 



Cl 



(fc) 
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^ci - X - \/3k/2. 



(12) 



writing the drive detuning as^ — lo^ — uJc- The dip in the 
heterodyne measurement of Fig. [I] appears within the bi- 
furcation region (Fig. [2k), indicating that this dip is the 
result of ensemble-averaging of the coherent heterodyne 
amplitude in the region of classical bistability. In Fig. ^p 
we see that the semiclassical and quantum simulation 
yield the same response outside the region of bistability. 
Within the region of bistability, quantum noise causes 
switching [36] between the two semiclassical solutions, 
one dim and one bright, with almost opposite phases. 
Therefore, averaging over the amplitudes of an ensemble 
of independent realizations gives rise to a destructive in- 
terference effect at the frequency where the amplitudes 
weighted by the switching rates are similar. Other forms 
of single atom bistability are known: single atom absorp- 
tive bistability |37j exists in the weak coupling regime 
in the good cavity limit, very different from the present 
discussion; more closely related is the single atom phase 
bistability of spontaneous dressed state polarization [31] 
which concerns the case where the atom and the cavity 
are in resonance 5 = 0, unlike our situation where the 
detuning 5 is the largest frequency scale apart from the 
cavity frequency. 

As the drive increases, and unlike the Kerr oscilla- 
tor, the frequency extent of the bistable region shrinks 
and eventually vanishes at the upper critical amplitude 
ic2 = g/V^- In the effective theory the upper crit- 
ical point C2 is located very close to the bare cavity 
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Figure 2. (Color online). Solution to semiclassical equa- 
tion (fTo|, using the same parameters as Fig. [lb. (a) Ampli- 
tude response as a function of drive frequency and amplitude, 
for same parameters as in Fig. [lb. The region of bifurcation 
is indicated by the shaded area, and has corners at the crit- 
ical points Cl, C2. The dashed lines indicate the boundaries 
of the bistable region for a Kerr oscillator (Duffing oscilla- 
tor) , constructed by making the power-series expansion of the 
Hamiltonian to second order in N/Nait- The Kerr bistability 
region ^5) matches the JC region in the vicinity of Ci but 
does not exhibit a second critical point, (b) Cut through (a) 
for a drive of 6.3^i, showing the frequency dependence of the 
classical solutions (solid blue). For comparison, the response 
for from the full quantum simulation of Fig. [lb is also plotted 
(dashed red) for the same parameters, (c) Cut through (a) for 
driving at the bare cavity frequency, showing the large gain 
available close to C2 (the 'step'). Faint lines indicate linear 
response, (d) Same as (c), on a linear scale. 
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Figure 3. (Color online) Symmetry breaking. State- 
dependent transition frequency versus excitation number: 
(a) for the JC model, parameters as in Figs. IT] and [2] (b) for 
the model extended to 2 qubits, 5i = —1.0 GHz, ^2 = 
—2.0 GHz, gi = g2 = 0.25 GHz; (c) for the model extended to 
one transmon qubit, tuned below the cavity, Ec ~ 0.2 GHz, 
Ej = 30 GHz, g = 0.29 GHz. In all panels, the transition fre- 
quency asymptotically returns to the bare cavity frequency. In 
(a) the frequencies within the ct^ = ±1 manifolds are (nearly) 
symmetric with respect to the bare cavity frequency. For (b), 
if the state of one ('spectator') qubit is held constant, then 
the frequencies are asymmetric with respect to flipping the 
other ('active') qubit. In (c), the symmetry is also broken 
due the existence of higher levels in the weakly anharmonic 
transmon. 



frequency. This indicates that for driving at the bare 
cavity frequency, there is no bistability, but rather a fi- 
nite region (a 'step') in the vicinity of the critical point 
(Fig. [2];) , where the response becomes strongly sensitive 
to the drive amplitude. The size of the step can be shown 
to be a factor of r = ^bright /^dim = 2g^ /k5 in amplitude, 
and represents a very high gain (dA/d^ = \f2gl iC'I'^b^l'^') 
in the strong-dispersive regime. Above the step we see 
that the response approaches the linear response of the 
bare cavity as A^ ~ Abarc- 



From the semiclassical Eqs. ( 10 ), (111 it follows that for 



A 3> 1 the response of the system will have an approx- 
imate symmetry of reflection with respect to the bare 
cavity frequency v4(S1,ct^ ~ +1) ~ A(— n, tj^ = —1). 
Therefore the response at the bare cavity frequency will 
be nearly independent of the state of the qubit. In or- 
der to translate the high gain available at the step into 
a qubit readout, it is necessary to break the symmetry 
of the response of the system between the qubit ground 
and excited states, such that the upper critical power ^ci 
will be qubit state dependent. In the JC model the sym- 



metry follows from the weak dependence of the decou- 
pled Hamiltonian H on the qubit state for high photon 
occupation. However, the experimentally-observed state 
dependence may be explained by a symmetry breaking 
caused by the higher levels of the weakly anharmonic 
transmon, or by the presence of more than one qubit, 
see Fig. [3] When designing a readout scheme that em- 
ploys such a diminishing anharmonicity, the contrast of 
the readout is a product of both the symmetry break- 
ing and the characteristic nonlinear response of the sys- 
tem near the critical point Ci- Experiments |14| were 
able to use this operating point to provide a scheme for 
qubit readout, which is attractive both because of the 
high fidelities achieved (approaching 90%, significantly 
better than is typical for linear dispersive readout in cir- 
cuit QED [351 [22]) and because it does not require any 
auxiliary circuit elements in addition to the cavity and 
the qubit. During preparation of this manuscript we be- 
came aware of a theoretical modeling of the high-fidelity 
readout by Boissonneault et al. |4Dj. 
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